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ABSTRACT 

We report on observations of the blazar W Comae (ONh-231) with the Solar Tower Atmospheric Cherenkov 
Effect Experiment (STACEE), a wavefront-sampling atmospheric Cherenkov telescope, in the spring of 2003. 
In a data set comprising 10.5 hours of ON-source observing time, we detect no significant emission from W 
Comae. We discuss the implications of our results in the context of the composition of the relativistic jet in W 
Comae, examining both leptonic and hadronic models for the jet. We derive 95% confidence level upper limits 
on the flux at the level of 1.5-3.5 x 10"'** cm~^ s"' above 100 GeV for the leptonic models, or 0.5-1.1 x 10"'° 
cm"^ s"' above 150 GeV for the hadronic models. 
Subject headings: gamma-rays: observations — BL Lacertae objects: 
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1. INTRODUCTION 

The current catalog of extragalactic TeV gamma-ray 
sources consists of blazars, which are among the brightest and 
most rapidly variable objects in the sky. The unusual observa- 
tional properties of blazars are usually explained in terms of 
relativistic bulk moti on of the e mitting region in a jet parallel 
to the line of sight iUrrv & Padovani. 19953 . The continuum 
emission of blazars is typically nonthermal and contains two 
broad peaks, one at lower energies (radio to X-ray) and one 
at higher energies (keV to TeV). Although the low-energy 
peak is generally assumed to result from synchrotron radia- 
tion from high-energy electrons in the jet, the mechanism for 
the production of the high-energy radiation is still a subject of 
debate, and several competing models exist. 

In "leptonic" jet models, the high-energy radiation is pro- 
duced by inverse Compton scattering of low-energy photons 
from a population of electrons and/or positrons; these mod- 
els are favored for the well - studied blazars Mrk 42 1 and 501 
dConni & Aharonianiri999t iK^^onelko et al J 12001. In the 
synchrotron self-Compton model (SSC), a single population 
of electrons and positrons produce the target photon field via 
synchrotron radiation, and subsequently upscatter the syn- 
chrotron photons to gamma-ray energies (Bloom & Marscher 
There may also be an external Compton (EC) 
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component of target photons from the accretion disk 
or ambient medium ( Dermer, Schlickeiser, & Mastichiadi^ 
1992; Sikora. Begelman. & Rees 1994). Alternatively, in 
"hadronic" jet models, protons play a central role. In hadronic 
models the high -energy radiation is attributed to photome- 
son processes (Mannheim '1993') or synchrotron radiation 
from protons or muons, as in the synchrotron proton blazar 
(SPB) mod el ( Muck e & Profli eroe 2000; Aharonian 200^ 
iMiicke etai.^200 3). The study of hadronic models was origi- 
nally motivated by the hypothesis th at the highest-ener gy cos- 
mic rays are produced in blazar jets fMannheim' 1995'). 

The BL Lac object W Comae (ON-i-231) may pro- 
vide an excellent test case for h adronic jet models 
JBoettcher. Mukheriee & Reimedl2()02l) . The transition be- 
tween the low-energy and high-energy peaks in the contin- 
uum of W Comae appears clearly in X-ray data taken by 
the BeppoSAX satelHte (Tagliaferri et al. 2000). These high- 
quality observations of the transition region place tight con- 
straints on leptonic models, requiring the predicted gamma- 
ray emission to cut off sharply above 100 GeV. In contrast, 
hadronic models may allow for significant emission above 
100 GeV. Observations by the EGRET detector aboard the 
Compton Gamma Ray Observatory show a hard power law 
spectrum (photon spectral index a= 1.73 ±0.18) exte nding 
up to about 10 GeV with no sign of any cutoff ( Hartma n et alJ 
Il999h . Yet the object has not been detected at energies above 
300 GeV, despit e repeated obser vation by the Whipple 10- 
m instrument (Horan et al. "2003"). At a redshift of 0.1, ab- 
sorption of gamma rays in this energy range by pair produc- 
tion 77 — > e'^e~ against the extragalactic background light 
(EBL) may be significant, but only at energies above about 
500 Ge V (Primack et al. 1999, 2001; Malkan & Stecker 2001] 
'Aharonian 2001 ), so the intrinsic emission spectrum of W Co- 
mae should be directly observable at energies lower than this. 

The Solar Tower Atmospheric Cherenkov Effect Experi- 
ment (STACEE) currently operates above an energy threshold 
of about 100 GeV for gamma rays. W Comae has been ob- 
served in the past with previous versions of the STACEE de- 
tector (Theoret 2001). New observations of W Comae were 
carried out in the spring of 2003 by the most recent version of 
the detector, STACEE-64. 
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2. THE STACEE DETECTOR 

STACEE is an atmospheric Cherenkov telescope which 
uses as its primary optic an array of solar mirrors (heliostats) 
at the National Solar Thermal Test FaciUty (NSTTF) in Al- 
buquerque, NM. STACEE uses the independently steerable 
heliostats to collect Cherenkov light from extensive air show- 
ers initiated by astrophysical gamma rays. Secondary mirrors 
are used to image the heliostat field onto a camera of pho- 
tomultiplier tubes (PMTs), producing a one-to-one mapping 
between heliostats and PMTs. Optical concentrators based 
on the DTIRC design dNing. Winston. & O'CiallagheJI 19871) 
widen the aperture of each PMT from 5 cm to 1 1 cm and re- 
strict its field of view to reduce the number of night sky back- 
ground (NSB) photons detected. High-speed electronics mea- 
sure the charge and relative aiTival times of the PMT pulses, 
and impose a coincidence trigger. The STACEE detector is 
thus a wavefwnt- sampling detector, measuring the intensity 
and arrival time of the naiTow wavefront of Cherenkov light 
at different spots on the ground. These measurements enable 
the offline reconstruction of the primary energy and shower 
aiTival direction. 

STACEE achieves a low energy threshold (~ 100 GeV) for 
the detection of gamma rays (compared with contemporary 
single-dish imaging Cherenkov telescopes) primarily because 
of the extremely large available mirror area (A ~ 10^ m^); 
the energy threshold scales approximately as A"'/^ ("Weekesj 
[l988). Other instruments using a similar solar tower design 
include the Cherenkov Low Energy Sampling and Timing 
Experiment (CELESTE) in Themis, France ( de N aurois et aL 
1^02), and the Solar Two detector in Barstow, CA jTripathil 
|2002). These are the only experiments currently operating in 
the Northern Hemisphere which have their peak sensitivity to 
gamma rays below 300 GeV. 

STACEE has undergone a series of upgrades and improve- 
ments over its lifetime. The prototype experiment, called 
STACEE-32, used 32 heliostats with a total mirror area 
of 1200 m^ (Hannaetal. 2002), and successfully detected 
gamma-ray emission from the Crab Nebula at a peak en- 
ergy of 190 GeV (Oseretal. 2001). An upgraded detector 
called STACEE-48 used an additional 16 heliostats (for a to- 
tal of 48) and employed improved high-speed programmable 
delay and trigger electronics (Martin & Ragan 2000); it de- 
tected gamma-ray emission from Mrk 421 during a period 
of intense flaring activity (Boone etal. 2002), as part of a 
multi-wavelength variability study. The final stage of con- 
struction, called STACEE-64, was completed in the summer 
of 2001. An additional 16 heliostats were instrumented, for a 
total of 64 heliostats with over 2300 m^ of miiTor area, filling 
in the geographical front and center of the heliostat field. The 
charge-integrating ADCs used by prior versions of STACEE 
were also removed, and replaced with waveform digitizers or 
"flash ADCs" (FADCs). A brief overview of these i mprove- 
ments follows; additional detail is available elsewhere dScalzd 
12004) . 

2.1. Optics upgrades 

Figure ^ shows the NSTTF heliostat field, with the 
STACEE heliostats indicated. The heliostats are chosen to 
provide the most uniform ground coverage possible, subject 
to constraints imposed by the crowding of heliostat images 
(and PMT apertures) in the focal plane of each secondary 
mirror In addition to the three secondaries (1.9 m diameter) 
used by STACEE-48, two additional secondaries (1.1 m di- 
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Fig . 1 . — Map of the heliostat field used by the STACEE experiment, show- 
ing updates to the experiment over time. Heliostats are numbered according 
to the trigger subcluster to which they belong. Clusters 0, 1,2, and 3 con'e- 
spond to the "east" and "west" cameras; most of these heliostats were used 
in the STACEE-32 prototype. Clusters 4 and 5 correspond to the "north" 
camera which was added for STACEE-48. Clusters 6 and 7 correspond to 
the new STACEE-64 "southeast" and "southwest" cameras. The base of the 
NSTTF solar tower, which houses the STACEE optics and electronics, marks 
the origin of the coordinate system. 



ameter) were built for STACEE-64 to allow use of additional 
heliostats at the front of the field. A set of 16 additional cam- 
era elements (PMT -i- DTIRC -i- housing) were instrumented 
for the new heliostats. All 64 camera elements were carefully 
calibrated in the lab and redistributed among the cameras ac- 
cording to measured differences in quantum efficiencies, in 
order to equalize the optical throughput on all detector chan- 
nels. 

Improved understanding of the heliostat optics via Monte 
Carlo simulations and measurements of the heliostat point- 
spread function using a CCD camera prompted adjustments 
to the optical elements. Each STACEE heliostat has 25 mir- 
ror facets which c an be independently focused and aligned 
('Hanna et al. 2002); the point spread function is periodically 
evaluated and adjustments made, if necessary, using a laser 
look-back method. Prior to the W Comae observing cam- 
paign, the facets on many heliostats were re-adjusted to fix 
the overall optical axis within the heliostat body frame to a 
standard position, thereby improving pointing accuracy. The 
absolute pointing of each heliostat was calibrated to an accu- 
racy of 0.05° using drift scans of bright stars. 

2.2. Trigger 

STACEE uses a two-lev el coincidence trigger described in 
detail in Martin & Raganl (.2000). PMT camera elements are 
grouped into subclusters of eight PMTs each. The signal from 
each PMT is AC-coupled, amplified, and discriminated before 
entering the subcluster trigger electronics. The coincidence 
trigger accounts for channel-to-channel delays among the dis- 
criminated PMT signals, which arise from the geometry of the 
shower wavefront and the differences in propagation times 
from different parts of the heliostat field; the programmed 
delays are accurate to 1 ns. A coincidence test among de- 
layed PMT signals is then made at the subcluster level, and 
a further coincidence among subclusters is required to trigger 
event readout. The number of coincident PMTs in a subclus- 
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ter, and the number of coincident subclusters, are chosen to 
optimize the quahty factor for the rejection of hadronic air 
showers according to Monte Carlo simulations. The discrim- 
inator threshold is then set at a level which makes the overall 
event trigger rate from chance coincidences of pulses due to 
fluctuations of NSB photons negligible (less than 0.2 min"'). 

The STACEE-64 trigger topology is shown in Figure^ For 
the W Comae observations, the discriminator threshold was 
set to 140 mV (about 5.5 photoelectrons); five out of eight 
PMTs in a subcluster were required to cross threshold within 
a narrow coincidence window to generate a subcluster trigger. 
The precise width of this coincidence window as applied to 
a given series of pulses varies between 8 and 24 ns due to 
the im plementation of the trigger logic (see Martin & Rag^ 
12000*): the mean width for a two-channel coincidence is 16 ns. 
Five out of eight subcluster triggers within a window of 16 ns 
were necessary to trigger event readout. 

2.3. FADCs 

The FADCs represent a major upgrade to the STACEE 
experiment. Access to the full digitized PMT waveform 
allows not only more accurate measurements of the tim- 
ing and intensity of the wavefront, but also measurement 
of charge-timing correlations, such as the distribution of 
Cherenkov photon arrival times at each heliostat. Vari- 
ous new methods which use FADC data to reject hadronic 
even ts are currently under study by the STA CEE collabora- 
tion (|Sciiioiitan|200l |Zwierm^^ The FADCs 
are also routinely used to calibrate and monitor the gains of 
the STACEE PMTs using a custom-designed laser calibration 
system (Hanna & Mukheriee 2001). 

The FADCs used by STACEE-64 are a commercial sys- 
tem produced by Acqiris, Inc. The system comes in modu- 
lar pieces, with four channels per board. Up to six boards 
(24 channels) fit into a special crate, which runs the Linux 
operating system. A real-time Linux driver for the system 
was developed by our group. Each FADC channel samples 
at 1 GS/s with a dynamic range of 1 V during normal astro- 
nomical observations. A sampled PMT signal from an actual 
Cherenkov event is shown in Figure|2] The zero points of the 
FADC inputs are calibrated to a precision of 1 mV RMS, and 
the channel-to-channel gains of the system are equalized to 
within 0.5%. 

In the 2001-2002 observing season, 32 FADC channels 
were instrumented, each sampling the signal from two PMT 
channels chosen so that the incoming Cherenkov signals 
would be well-separated in time. By the 2002-2003 observ- 
ing season, 64 FADC channels were instrumented so that each 
FADC sampled the signal from exactly one PMT. 

3. OBSERVATIONS AND QUALITY CUTS 

STACEE-64 was used to observe W Comae in the spring 
of 2003, from late March to early June. To maximize sensi- 
tivity and robustness to systematic effects, ST ACEE operates 
in an "ON-OFF" integration mode (lOser et al . 2001). In ON- 
OFF observation, each gamma-ray source is assigned another 
region of the sky at the same declination which contains no 
known gamma-ray source. The two regions are usually sepa- 
rated by 30 minutes in right ascension. STACEE then tracks 
the OFF region immediately before or after the ON region, 
so that the same range of azimuth and elevation are observed, 
and the detector sensitivity (which depends on azimuth and 
elevation) is matched in both halves of the pair. Thirty-four 
ON-OFF pairs were taken on W Comae, for a total of 13.5 


















20 40 60 80 100 120 140 160 180 

Time (ns) 



Fig. 2. — Digitized PMT signal from a typical air shower event. Dot-dash 
line: FADC baseline. Dashed line: Discriminator threshold. Dotted line: 3 
standard deviations above baseline for NSB fluctuations. 



hours of ON-source observing time. 

Several quality cuts were then imposed on the data set, as 
described below. These quality cuts are to be understood as 
"pairwise time cuts", in the sense that if data in a certain time 
interval during one run are flagged and removed from the data 
set, then data in the corresponding time interval in the other 
run of the pair (i.e., spanning the same local horizon coordi- 
nates) are also removed. 

The control system for the NSTTF heliostats continually 
records status information which can be used to pinpoint mal- 
functions on particular heliostats. This status information is 
logged nightly and merged with the regular STACEE data 
product. Portions of runs during which any heliostat malfunc- 
tioned were flagged and removed from the data set. 

As observed in Boone et al. (2002), the subcluster trigger 
rates are very sensitive to fluctuations in sky conditions, since 
they vary as a high power of the discriminator rates on PMTs 
within the cluster, which in turn are dominated by fluctua- 
tions of NSB photons. Large deviations of the subcluster trig- 
ger rates from a smooth linear trend have been seen to cor- 
relate strongly with unstable weather conditions. After he- 
liostat cuts, an iterative procedure was used to identify and 
flag affected portions of the data. Each pair was divided into 
30-s time bins, and a least-squares fit to the time-averaged 
subcluster rates performed repeatedly, excluding outliers be- 
yond a certain tolerance at each iteration. In order for a bin to 
merit exclusion from the data set, deviations in the rates had 
to appear for at least four out of eight subclusters. After un- 
dergoing this process, if a pair of data had less than a minute 
of data remaining, this pair was simply excluded from further 
analysis. 

Additionally, some data were removed because of malfunc- 
tions in the data acquisition program resulting in the loss or 
corruption of FADC data, which is needed for the final anal- 
ysis (see below). About 2% of the original data set was re- 
moved for this reason. 

Finally, even though the pairwise time cuts equalize the ob- 
serving times in the ON-source and OFF-source data sets, a 
livetime correction is necessary. The time required for the 
data acquisition software to read out each triggered event ere- 
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ates dead time within the observing intervals, which is rate- 
dependent and measured electronically by the detector as the 
data are acquired. 

After all quality cuts, a total of 10.5 hours of ON-source ob- 
serving time remained, distributed among 32 ON-OFF pairs 
(see Table [Hi. Using the Li & Ma ( 1983) expression for the 
significance gives a (raw) positive significance of 4.6 standard 
deviations, at a livetime-corrected excess event rate of 4.8 ± 
1 . 1 min~' . 

4. SKY BRIGHTNESS CORRECTIONS 

Before deriving an integral flux or upper limit, it is very 
important to correct the raw trigger rate for any differences 
in NSB fluctuations between the ON-source and OFF-source 
fields. NSB photons entering the detector frequently pro- 
mote sub-threshold hadronic air showers above the detection 
threshold, even under the best conditions. The promotion rate 
increases as the sky brightness increases, so that if the ON- 
source field is slightly brighter than the OFF-source field, ei- 
ther due to stars or to subtle variations in atmospheric con- 
ditions, the resulting excess in the hadronic trigger rate can 
mimic a gamma-ray excess. Similarly, if the OFF-source field 
is brighter than the ON-source field, an existing gamma-ray 
excess may be masked; a net deficit may even appear An 
accurate measurement of the true gamma-ray rate therefore 
requires a correction for the promotion effect. 

We have developed three independent methods to correct 
for the effects of varying levels of NSB light on the coinci- 
dence trigger. We find that the excess we see from W Comae 
is consistent with the promotion of events representing only 
hadronic showers, with no additional gamma-ray component. 

4.1. Direct measurement of the hadronic promotion trend 

One method of correcting for field brightness effects was 
demonstrate d by the analysis o f the STACEE-48 observations 
of Mrk 421 ("Boone et al.''2002). The excess brightness in the 
STACEE ON-source field for Mrk 421 was due mainly to a 
single bright star, HD 95934 (magnitude 6.16 in the Z?-band). 
No FADC data were available at that time, but the promotion 
rate was estimated simply by taking ON-OFF pairs on another 
star, HIP 80460, with a similar Z?-band magnitude. The mea- 
sured promotion rate was then subtracted from the total excess 
to give the gamma-ray excess. 

Similar observations have been made for STACEE-64 and 
are depicted in Figure |3l Three stars of different magni- 
tudes and suitable declinations (21 Com, l CrB, and HIP 
89279) were selected and observed as if they were gamma-ray 
sources. The excess rate is plotted as a function of a weighted 
average current difference, 

where A/, is the ON-OFF difference in anode current on PMT 
j, Wj is the measured fraction of triggered events in which 
PMT j crossed threshold, and the angle brackets denote time- 
averaging over the entire data set. The promotion trend is well 
fit (x^/d.o.f. = 1.04) by a straight Une with slope 3.7 ±0.3 
min"' /iA"'; it represents a zero-signal baseline for gamma- 
ray sources observed by STACEE. 

The raw rate measurement for W Comae is also shown on 
the plot, and it is consistent with the promotion trend. Sub- 
tracting the expected number of promotion events from the W 
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Fig. 3. — Excess event trigger rates as a function of the characteristic cur- 
rent difference for sources observed in 2002-2003 by STACEE-64. Squares: 
star observations; circle: W Comae. Filled symbols show the raw excess, 
and open symbols show the excess after the dynamic threshold brightness 
correction. The lines show the best-fit linear trend constrained to go through 
the origin (solid) with 68% confidence errors on the slope (dashed), and the 
zero-rate baseline (dotted). 



Comae excess based on its characteristic current leaves a net 
excess of 0.2 ± 1 .2 min"' . 

The direct measurement of the promotion trend is intuitive 
and straightforward, and can be used to quickly assess how 
much of an observed excess is due to promotion. However, 
the promotion trend as measured reflects the observing sea- 
son on average, and in particular observing conditions dur- 
ing the star runs, rather than conditions which might pertain 
for individual runs in the data set for a particular gamma-ray 
source. An event-by-event, channel-by-channel treatment is 
therefore more desirable when deriving scientific results such 
as fluxes or upper limits. A sizable data set of star observa- 
tions proved quite valuable for testing such brightness correc- 
tion techniques. 

4.2. Software padding using waveform libraries 

Another technique, commonly referred to as padding, has 
been in use in some form by atmospheric Cherenkov tele- 
scopes since their inception. For example, in the early days of 
the Whipple 10-m imaging telescope, NSB conditions were 
equalized in hardware by using a small LED on the face of 
each PMT in a feedback loop (Weeke s et al.lll989h . Later 
the equalization was done in software, by adding random 
Gaussian deviates to the recorded ADC values to mimic fluc- 
tuations from increased NSB levels dCawlev 1993), and by 
then reimposing an additional selection criterion similar to the 
hardware trigger. For STACEE, the trigger depends not only 
on charge distributions but also on relative timing, since the 
light from an air shower reaches different parts of the detec- 
tor at different times. Knowledge of the PMT waveforms as 
provided by the FADCs is therefore necessary to do software 
padding. 

A software padding scheme has been implemented and 
tested by STACEE as follows: FADC waveforms correspond- 
ing to varying NSB levels were sampled from 16 of the PMTs 
in situ under controlled conditions. The waveforms were char- 
acterized by their RMS fluctuations and stored in a library 
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Fig. 4. — Excess event trigger rate as a function of FADC analysis tliresli- 
old in tlie library padding analysis of a subset of data taken on the star i CrB. 
The trigger condition is described in the text. Upright triangles: No padding 
of traces before retiiggering. Inverted triangles: Traces are padded to equal- 
ize RMS before retriggering. Dashed line: zero promotion excess. The figure 
shows that the minimum thi'eshold necessary to remove all hadronic promo- 
tion events from the library-padded data set is 150 mV. 



to be used in padding. In the padding analysis, these library 
waveforms were added to the waveforms in the observation 
data set so as to equalize the RMS fluctuations in both runs of 
a pair, as measured by the 400-ns section of the waveform im- 
mediately preceding the event trigger. An offline reimposition 
of the trigger required five out of eight PMTs in a subcluster, 
and five out of eight subclusters, to fire within a coincidence 
window of width 16 ns. 

When reimposing the trigger in software, the analysis 
threshold at which the FADC waveforms are discriminated 
must be increased over the nominal hardware discrimina- 
tion threshold. Adding background noise to waveforms from 
recorded events in a low-noise run equalizes only the proba- 
bility that an event will pass additional offline cuts relative to 
an event in the corresponding high-noise run. It is impossible 
to recover events which would have appeared in the low-noise 
data set if the additional noise had been physically present at 
the outset. A subset of the ON-OFF data for i CrB were used 
to optimize the offline analysis threshold (see Figure]^. The 
analysis threshold was set to 150 mV for all FADC channels 
based on the results of this optimization. 

The results of the library padding brightness correction pro- 
cedure are shown in Table For all three stars and for W 
Comae, the ON-OFF difference drops from a significant ex- 
cess in the uncorrected data to an insignificant difference in 
the corrected data. In particular, the significance drops from 
4.6 to 0.9 standard deviations for W Comae. 

4.3. Dynamic thresholds 

Although software padding has clearly been successful in 
eliminating the promotion excess, it has the drawback that 
the software discriminator threshold used to analyze FADC 
waveforms, and hence the energy threshold for the analysis, 
must be raised significantly. The desire to maintain a low 
discriminator threshold in the offline analysis of FADC data 
formed the motivation for investigating what we call the dy- 
namic thresholds technique. 

In the absence of NSB fluctuations, the distribution of 
pulse heights for pulses within the trigger coincidence win- 
dow on each channel would depend only on the intensity 
of Cherenkov light from air showers. The presence of NSB 
fluctuations on a channel has two effects on the pulse height 



spectrum. First, the linear superposition of fluctuating back- 
ground traces upon the underlying pulses from Cherenkov 
photons produces an additional statistical uncertainty in the 
pulse height measurement for each event, by changing the 
pedestal from which the pulse height is measured. NSB fluc- 
tuations therefore smear out the underlying air shower pulse 
height distribution. Second, NSB fluctuations occasionally 
cross the discriminator threshold within the coincidence win- 
dow even when no signal pulse is present, resulting in a sec- 
ond component to the pulse height distribution. 

The dynamic thresholds technique models both of these ef- 
fects quantitatively to predict changes in the analysis thresh- 
olds. When comparing the ON-source and OFF-source runs 
of a pair, the run with larger fluctuations in the baseline trace 
will have a higher promotion rate, and therefore a lower ef- 
fective discriminator threshold. To compensate for promo- 
tion, for each channel the threshold in the noisier run is raised 
so that the rate of background pulses (NSB plus hadronic 
Cherenkov light) within the coincidence window, and with 
pulse height above the analysis threshold, is the same in both 
halves of the pair When a coincidence condition among 
pulses is imposed offline, this prescription should remove a 
number of events equal to the expected number of promoted 
hadronic shower events, while retaining sensitivity to any pos- 
sible gamma-ray excess. 

For the OFF-source run, only hadronic air showers and 
NSB photons contribute to the pulse height distribution, so 
that the rate of pulses above the analysis threshold can be di- 
rectly measured. The pulse height distribution for the ON- 
source run, however, may also contain Cherenkov pulses from 
gamma-ray air showers. Thus, the form of the ON-source 
pulse height spectrum without the gamma-ray contribution 
must be predicted from the measured OFF-source pulse height 
spectrum and measurements of the NSB fluctuations. 

We used a semi-analytic approach to predict the form of 
the O N-source pulse height spectrum (for details, see Scalz^ 
\m)4t . Histo grams of individual waveform samples were ac- 
cumulated using the part of each waveform immediately pre- 
ceding the event trigger These histograms characterized the 
distribution of fluctuations, denoted .son(V) and .soff(^), in 
the waveform baseline (and hence in the effective discrimina- 
tor threshold) due to NSB fluctuations in each run. The pulse 
height distribution for all pulses within the coincidence trig- 
ger window for the OFF-source run, denoted ^off(V), was 
also measured. We predicted the form of the ON-source pulse 
height spectrum /?on(V) as follows: 

1. Fit .soff(V), ■son(V) and /?off(V) with smooth analytic 
forms to interpolate between thresholds. 

2. Normalize .soff(V) with respect to /?off(V) using the 
measured OFF-source PMT rate (which is dominated 
by NSB fluctuations) and subtract it from /?off(V) to 
obtain a pulse height spectrum representing only pulses 
due to hadronic air showers. 

3. Deconvolve this hadronic pulse height spectrum by 
■soff(V) to obtain ^o(V), the expected hadronic air 
shower pulse height spectrum in the absence of NSB 
fluctuations. 

4. Convolve Rq{V) by son(V) to account for the NSB fluc- 
tuations in the ON-source run. 

5. Normalize ioN(^) using the measured PMT rate ON- 
source and add it to the convolution of Rq^V) and 
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Table 1. Results of held brightness corrections 



50 100 150 200 250 300 350 400 450 500 

Pulse Height (mV) 

Fig. 5. — Pulse height distiibution of tiiggered events from a 10-minute 
segment of a typical OFF-source run on the star l CrB. The best-fit form of a 
Poisson distribution (expected NSB) plus an exponential is shown (solid line), 
as well as the predicted form for the ON-source run (dashed line) resulting 
from additional NSB fluctuations. The ON-source prediction is normalized to 
the same detector livetime as the OFF-source best fit. Horizontal dashed line: 
number of OFF-source pulses per mV at the nominal hardware discriminator 
threshold. Vertical dashed lines: analysis thresholds ON and OFF necessary 
to equalize the number of pulses above threshold over the course of the run. 



■son(V) to obtain the predicted form of /?on(V). 

Figure|5lshows an example of the results of this procedure for 
one STACEE channel in the data set for the star t CrB. A vi- 
sual inspection of Figure|5lindicates that the threshold on this 
channel must be increased by 15 mV. This figure represents 
the worst-case scenario, corresponding to conditions on one 
of the noisiest channels during observations of the brightest 
of the three padding stars. For the W Comae observations, 
the mean increase in the software discriminator threshold is 
approximately 2 mV, considerably smaller than the 10 mV 
adjustment on all channels necessary for software padding. 

The dynamic thresholds for a run were always calculated 
after quality cuts to eliminate systematics due to rapidly fluc- 
tuating sky conditions. To improve robustness against slow 
drifts in sky conditions, each 28-minute run was partitioned 
into three equal segments 560 seconds in length. The statis- 
tics within each segment were sufficient to establish the form 
of the sampled distributions, except in cases where time cuts 
had already removed most of the segment. If statistics within 
a segment did not permit a reliable prediction of the threshold 
adjustments, the segment and its counterpart in the other half 
of the pair were removed from the data set. Less than 3% of 
the quality-cut data set was removed in this way; the remain- 
ing data set comprised 10.3 hours of ON-source observing. 

The results of the dynamic threshold analysis for W Comae, 
and for the three stars observed by STACEE for promotion 
studies, are shown in Tabled As with library padding, no 
statistically significant excess remains on any of the stars, or 
on W Comae, after the dynamic threshold technique is applied 
to the raw data sets. 

5. FLUX LIMIT DETERMINATION 

To calculate an upper limit on the flux of gamma rays 
from W Comae, two main ingredients are necessary: an as- 
sumption about the shape of the spectrum, and knowledge of 
the detector response. The latter was produced from Monte 
Carlo simulations. We use the CORSIKA simulation package 
(|Heck et al. 1999) to model both hadronic and gamma-ray air 
showers. The detector was modeled using software written 
by the collaboration, including a full ray-trace of the detec- 
tor optics and detailed simulations of the discriminators, the 



Source 


tON'' (S) 


tOFF'-' (s) 


NON 


NOFF 


Sig. {at 


QuaUty cuts only: 












iCrB 


11777.4 


12064.9 


75505 


73060 


-1-10.99 


21 Com 


6053.2 


6148.7 


32354 


31299 


-1-6.16 


HIP 89729 


19087.2 


19191.8 


74823 


73292 


-1-5.03 


W Comae 


37161.0 


37955.5 


219031 


217085 


-1-4.60 


Library padding: 












L CrB 


mil A 


12064.9 


52423 


53367 


-Hi. 02 


21 Com 


6053.2 


6148.7 


22859 


23269 


-0.23 


HIP 89729 


19087.2 


19191.8 


53463 


53678 


-1-0.24 


W Comae 


37767.0 


37955.5 


155722 


156006 


-1-0.88 


Dynamic thresholds: 












t CrB 


11234.1 


11508.5 


58566 


60556 


-1.60 


21 Com 


5522.8 


5613.9 


25172 


25305 


-1-1.25 


HIP 89729 


18731.9 


18835.3 


61990 


62064 


-1-0.76 


W Comae 


36982.9 


37168.9 


182433 


182942 


-hO.67 



Note. — Results of field brightness connections for three stars of varying 
magnitude, and for W Comae. The live times and raw excesses of events in 
the ON and OFF runs, and the corresponding significance, are shown for each 
technique. 

"Integrated live time after quality cuts and dead time correction. Note that 
these represent equal amounts of actual observing time in the ON and OFF 
runs, but that differences in the trigger rate may lead to differences in the 
detector dead time between the ON and OFF runs. 

^Simificance in standard deviations (tr). Calculated according to lLi & m3 



delay and trigger systems, and the FADCs. The full simula- 
tion chain successfully reproduces important detector-related 
quantities, such as the discriminator rate on each PMT, the 
trigger rate of cosmic-ray air shower triggers at zenith, and 
the increase in the hadronic trigger rate resulting from NSB 
fluctuations. 

The response of the STACEE detector is characterized by 
the effective area, defined as a function of energy by 



A(E)= / P{x,y,E)dxdy 



(2) 



where P{x,y,E) is the probability that a gamma-ray air shower 
with shower core landing at position ix,y) produces an event 
trigger, and G is a region on the ground (the xy-plane) con- 
taining all (x,y) for which P{x,y,E) j^Q. The functional form 
of the effective area used to obtain the integral flux limits in 
this paper is an average of the results of simulations done in 
different regions of the sky, at a declination of -1-28.23° spaced 
in 10° increments in hour angle, each weighted by its expo- 
sure in the 2003 data set. The sensitivity to very low-energy 
air showers is highest for a source at transit. However, about 
80% of the 2003 W Comae data were taken between hour an- 
gles of H-20° and H-40°(zenith angles between 10° and 26°). 
The large zenith angles result in a loss of sensitivity for low- 
energy air showers with respect to that expected for the source 
at transit. First, the Cherenkov photon density of the shower 
decreases, and the area of the Cherenkov light pool increases, 
with increasing zenith angle. Second, off-axis aberrations in 
the heliostat optics broaden the heliostat point spread func- 
tion for sky locations far from transit, decreasing the overall 
optical throughput. 

The effective area must also account for any energy depen- 
dence introduced by the offline trigger cut used in the anal- 
ysis. Of the three analysis methods available to correct for 
field brightness differences, the dynamic threshold technique 
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Fig. 6. — STACEE-64 simulated trigger rate vs. energy for a power law 
differential flux spectrum of gamma rays with photon index a = 1.55 (in- 
verted triangles), a = 1.73 (circles), and a = 1.91 (triangles). All spectra are 
normalized to the same integral flux above 50 GeV. 



maintains the best sensitivity at low energies, and we use this 
technique to estabHsh our upper Hmit. 

In previ ous S TACEE pubHcations ("Oser et aTl 1200 It 
iBoone et al]|2nn3), the energy threshold has been quoted as 
the energy at which the differential gamma-ray trigger rate 
(effective area times expected photon flux from the assumed 
model) reaches a maximum. We continue to do so in this pa- 
per, quoting an energy threshold for each assumed model ac- 
cordingly. We assume that the errors in the absolute through- 
put calibration of the independent elements (heliostat optics, 
PMT quantum efficiency, etc.) of the STACEE detector are 
Gaussian and uncorrected. These errors may then be added 
in quadrature to obtain the systematic uncertainty in the ab- 
solute energy scale calibration, which we estimate to be 20%. 
The systematic errors on the energy threshold are determined 
by rescaling the energy scale at which the effective area is 
evaluated by ±20%, and calculating the resulting shift in the 
maximum of the differential gamma-ray trigger rate. 

Figure|6ldemonstrates the calculation of the energy thresh- 
old for an extrapolation of the EGRET power law spectrum 
for W Comae {a = 1.73 ±0.18), using the exposure-averaged 
effective area. Curves representing the low and high values 
for the spectral index within the EGRET experimental errors 
(a = 1.55 or 1.91) are also shown, normalized to the same 
integral flux above 50 GeV. At transit, at a zenith angle of 
6.73°, the energy threshold is about 120 GeV for a = 1.73. 
The energy threshold using the exposure-averaged effective 
area curve is 170 GeV. For the a= 1.73 model, therefore, we 
quote our result at an energy threshold of 170 ±40 GeV. This 
threshold changes by about 20 GeV if the low or high value 
of a (1.55 or 1.91) is used instead. 

Tables HE] and g] show upper limits from STACEE ob- 
servations for selected model s of the g amma ra y emis sion 
from W Comae dBo ettcher. Mukheriee & Reimen l2002h . In 
each case the optical depth for pair production against the 
EBL has been taken into account. The models are based 
on the semi-an alytic approaches of Pri mack et alJ ( 1999) and 
iPrimack et alJ 12001 ). which use a top-down, hierarchical 
treatment of structure formation as tuned to fit recent EBL 
measurements. The qualifiers "low" and "high" refer to the 
extreme low (LCDM-Salpeter; Primack et al. 1999) and high 
(Kennicutt; Primack et alJl200 1 ) predictions for the optical 
depth to pair production at z = 0. 102 from s pecific instances of 
these semi-analytic models, as detailed in lAharonianl j2001l) 



Table 2. Integral flux upper limits for power law spectra 



Emission Model r^-^" f'thrcsh'' 95% CL'^ *(> f'thrcsh)'' 



a =1.55 low 191!|^ 0.48 7.28 

a =1.55 high 186!^ 0.52 3.87 

a =1.73 low 186!^J 0.49 2.23 

a =1.73 high HO^fS 0.60 1.49 

a = 1.91 low 168!38 0.56 0.73 

a= 1.91 high 160!^ 065 0.52 



Note. — Upper limits on the integral flux of gamma rays from W Co- 
mae, assuming that the differential flux of photons follows a power law 
(^( E) oc E^°'). Power la ws were extrapolated from the Third EGRET Cata- 
log I Ha rtman et aH199^ and corrected for EBL absorption. The mean spec- 
tral index o = 1.73 is used, as well as low and high values consistent with the 
eiTors on the power-law fit to the EGRET data. 

"Optical depth to pair production 77 — > e+e" against the EBL from a red- 
shift of z = 0. 102, from Aharoniat) 120011) . The extreme low and high predic- 
tions for r.^-y were used in this table. 

''STACEE energy threshold in GeV for the assumed model. The energy 
threshold is defined as the energy at which the differential trigger rate (ef- 
fective area times gamma-ray flux) reaches a maximum. Systematic errors 
as shown represent propagation of the error in the calibration of the absolute 
energy scale using the mean effective area curve. 

'^STACEE 95% confidence level upper limit on the integral flux above the 
energy threshold, in units of 10"'" cm"^ s"'. 

''Expected flux above the energy threshold according to the emission model, 
in units of 10"'" cm"^ s"'. 



Table 3. Integral flux upper limits for leptonic models 



Emission Model 


^thresh 


95% CL 


*(> ^thresh) 


SSC fit 1 


1 1 3+25 


1.35 


0.0278 


SSC fit 2 


1.39 


0.0603 


SSC fit 3 


108!?^ 


1.55 


0.0876 


SSC fit 4 


07+24 


2.12 


0.139 


SSC fit 5 


86+'5 


3.49 


0.17 


SSC fit 6 
SSC fit 7 


95+°3 
100+1 

112+26 


2.36 
1.94 


0.0276 
0.331 


SSC+EC fit 8 


1.37 


0.0303 


SSC+EC fit 9 


110 J 


1.41 


0.0344 


SSC+EC fit 10 


108:^7 


1.44 


0.0467 



Note. — Upper limits on the integral flux of gamma rays from W Comae, 
according to various leptonic (SSC, SSC+EC) models of the relativistic jet 
described in Boettcher. Mukheriee & Reimer 1 2002), denoted BMR02 above. 
Due to the sharp cutoff in the expected gamma-ray emission spectra above 
~ 100 GeV, EBL absorption should be negligible for the leptonic models. 
The table column headings are defined in the notes to Tablel2l 



Table 4. Integral flux upper limits for hadronic models 



Emission Model 




^thresh 


95% CL 


*(> ^thresh) 


SPB fit 1 


low 


,97+22 

-\i 


0.93 


0.87 


SPB fit 1 


high 


126!j| 


1.05 


0.66 


SPB fit 2 


low 


150^1^ 


0.72 


2.48 


SPB fit 2 


high 




0.77 


1.91 


SPB fit 3 


low 


164i 


0.58 


0.40 


SPB fit 3 


high 


157+22 

'-"-26 


0.67 


0.26 


SPB fit 4 


low 


163+30 


0.59 


0.44 


SPB fit 4 


high 


0.69 


0.30 


SPB fit 5 


low 


160^" 


0.60 


0.46 


SPB fit 5 


high 


141+36 

'^'-23 


0.82 


0.36 



Note. — Upper limits on the integral flux of gamma rays from W Comae, 
according to various hadronic (SPB) models of the relativistic jet described 
in Boettcher, Mukheriee & Reimer ( 2002), denoted BMR02 above. The table 
column headings are defined in the notes to Tablel2l 
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E (GeV) 

Fig. 7. — Integral flux in cm"^ s"' for each of the five SPB models from 
|Boettcher. Mukheriee & Reimei 1 2002 ), and for the extrapolation of the best- 
fit EGRET power law ( Hartman et al. 1999). Model predictions for both low 
and high EBL optical depth are shown, along with STACEE upper limits for 
mean EBL optical depth. Upper limits fr om observations w ith the Whipple 
10-meter from 1995-2000 are also shown IHoran et al j200^ . 



and Boettcher. Mukheriee &Reimer(2002). 

It is apparent that, even when EBL absorption is taken into 
account, the STACEE resuhs place strong constraints on the 
extrapolation of the EGRET power law to higher energies. A 
steeper value of the spectral index, within the uncertainties, is 
favored. As noted by Boettcher, Mukheriee & Reimer ( 2002), 
however, the EGRET data were not taken simultaneously with 
the measurements at other wavelengths. The possible power 
law spectrum may therefore represent a transient flaring state. 
Moreover, the EGRET power law is fit to a co-addition of 
many low significance (~ 2 standard deviations) detections 
from different epochs, with large uncertainties on the flux and 
spectral index. In fact, the leptonic model fits shown used 
only the Beppo-SAX X-ray data, which is of higher quality; 
no pure SSC model can fit both the X-ray and EGRET data 
simultaneously. 

The STACEE observations are unable to place interest- 
ing constraints on most of the leptonic models, with the ex- 
pected fluxes being an order of magnitude or more below the 
STACEE limit in each case. The low predicted fluxes may be 
anticipated from the sharp spectral cutoff in each model at or 
below 100 GeV, due to a similar sharp cutoff in the electron 
energy spectrum at 7c^niax 5 x 10"*. Since the dominating 
factor is the availability of energy in relativistic electrons, the 
cutoff energy in the modeled gamma-ray spectra is approxi- 
mately independent of the density of external target photons in 
SSCh-EC models. The cutoff in the electron energy spectrum 
is constrained tightly by the X-ray data, which lie between the 
falling edge of the low-energy (electron synchrotron) peak in 
the spectrum and the high-energy peak. 

On the other hand, most of the hadronic models predict inte- 
gral fluxes above the energy threshold which are only slightly 
below the corresponding upper limits from the STACEE ob- 
servations. In particular, we can exclude SPB model 2 at 95% 
confidence. Among all the hadronic models shown, model 



2 has the highest cutoff for the proton energy distribution 
(7p,max = 3 X 10^), yet the magnetic field, the Doppler factor 
of the jet, etc., are comparable to the other models. 

Figure shows the STACEE upper limits for the hadronic 
models in a graphical form. The figure shows that STACEE is 
quite close to the sensitivity necessary to begin to distinguish 
between the various hadronic models for W Comae. How- 
ever, the constraints provided by precise, simultaneous ob- 
servations at other wavelengths are necessary, in addition to 
STACEE observations, because of the potential high-energy 
variability of blazars such as W Comae. The constraints 
placed by the X-ray observations from SAX were instrumen- 
tal in shaping the predictions for the gamma-ray spectra for 
the hadronic models presented here. Figure Q also shows 
upper limits from high-energy observations made using the 
Whipple 10-meter imaging Cherenkov telescope (Horan et aQ 
|2003) from 1995 to 2000. (The energy thresholds for the 
Whipple limits were each calculated assuming a power law 
with spectral index a = 2.5, and do not yet address specific 
emission models.) Each instrument operates within a differ- 
ent regime of sensitivity and energy, and the scientific pay- 
off for a variable target is maximized when observations are 
made simultaneously. Further STACEE observations in the 
future, coupled with new offline hadronic rejection techniques 
to improve the flux sensitivity, should allow STACEE to fur- 
ther explore the parameter space for the hadronic models, es- 
pecially when contributing to a multi-wavelength monitoring 
campaign. 

6. CONCLUSIONS 

The STACEE-64 observations in the spring of 2003 were 
made at a lower energy threshold than any other atmospheric 
Cherenkov telescope has yet attained for W Comae. STACEE 
detects no significant emission from W Comae, resulting in 
95% CL upper limits on the integral flux above this threshold 
in various hadronic emission models at the level of 10"'" cm"^ 
s~'. While leptonic models predict a flux which falls below 
this level, extrapolations of the best-fit EGRET power law, 
and some synchrotron-proton hadronic models, predict an in- 
tegral gamma-ray flux above the energy threshold close to, or 
exceeding, the upper limit from STACEE observations. Ad- 
ditional STACEE observations planned for the spring of 2004 
ought either to exclude these models at a significantly higher 
confidence level, or to detect gamma-ray emission from W 
Comae if these models provide an adequate description of the 
source. 
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